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Abstract: An angular acceleration sensor can be used for the dynamic analysis of human 
and joint motions. In this paper, an angular acceleration sensor with novel structure based on 
the principle of electromagnetic induction is designed. The method involves the construction 
of a constant magnetic field by the excitation windings of sensor, and the cup-shaped rotor 
that cut the magnetic field. The output windings of the sensor generate an electromotive 
force, which is directly proportional to the angular acceleration through the electromagnetic 
coupling when the rotor has rotational angular acceleration. The mechanical structure and the 
magnetic working circuit of the sensor are described. The output properties and the 
mathematical model including the transfer function and state-space model of the sensor are 
established. The asymptotical stability of the sensor when it is working is verified by the 
Lyapunov Theorem. An angular acceleration calibration device based on the torsional 
pendulum principle is designed. The method involves the coaxial connection of the angular 
acceleration sensor, torsion pendulum and a high-precision angle sensor, and then an initial 
external force is applied to the torsion pendulum to produce a periodic damping angle 
oscillation. The angular acceleration sensor and the angle sensor will generate two 
corresponding electrical signals. The sensitivity coefficient of the angular acceleration 
sensor can be obtained after processing these two-channel signals. The experiment results 
show that the sensitivity coefficient of the sensor is about 17.29 mv/Krads . Finally, 
the errors existing in the practical applications of the sensor are discussed and the 
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corresponding improvement measures are proposed to provide effective technical support 
for the practical promotion of the novel sensor. 

Keywords: sensor; angular acceleration; electromagnetic induction; mathematical model; 
calibration; error analysis 



1. Introduction 

Angular acceleration is a common physical quantity in a rotating system, which can reflect the 
vibrational state of the rotation angle for the rotation shaft and transmission equipment while working. 
The responses of the rotation shaft and transmission equipment to various stimulations can be analyzed 
by measuring the angular acceleration, and especially the dynamic interference of rotary systems can 
be displayed in the form of angular acceleration. The measurement of rotating mechanical angular 
acceleration has wide applications, including the fields of automobiles, the military, aerospace, 
industry, electronics, etc. [1-5]. 

Angular acceleration can be measured indirectly by employing rotation angle sensors or angular 
speed sensors. However, the signal processing of these methods are very troublesome. Especially, the 
problems of delay characteristics and noise amplification are hard to resolve, which has prompted the 
development of direct angular acceleration measurement techniques. The applications of linear 
acceleration sensor are very wide [6-9], although it is more urgent to develop an angular acceleration 
sensor without the range limitation of the rotation angle. There exist the following problems in the 
presented research and development of angular acceleration sensors: lack of range limitation of the 
rotation angle, small phase lag, high signal-to-noise ratio and a wider band [10]. 

In recent years, some experts and scholars have conducted a series of researches on angular 
acceleration measurement. Moody [11] proposed a new sensitive angular accelerometer, which 
contained a superconducting test mass suspended by a weak flexure pivot; a superconducting quantum 
interference device amplifier and superconducting circuit were used to detect the angular displacement 
of the test mass. Sparks [12] described an electroformed micro-machining technology capable of 
producing both angular rate sensors and accelerometers, which were fabricated by CMOS and 
provided improved signal output. Jianli Li [13] proposed a novel micro electromechanical system 
(MEMS) pendulum angular accelerometer with electrostatic actuator feedback; it adopted a proof 
pendulum with optimized moment of inertia, suspended to dual anchors by a pair of torsion spring 
beams, as sensing component. Tomikawa [14] dealt with a piezoelectric angular acceleration sensor 
which can accurately distinguish a signal of rotational motion from other linear motion ones. In order 
to overcome the influence of gravity in angular acceleration measurements, a system for knee joint 
motions was constructed with a couple of accelerometers by Yamamoto [15]. The accurate angular 
acceleration could be calculated as the difference between two output values of the accelerometers. A 
bulk silicon micro-machined structure 6.4 mm x 6.4 mm in size designed and fabricated by Mizuno 
could detect angular acceleration and acceleration simultaneously and independently [16]. An angular 
acceleration measurement system without moving parts using a PVDF and a piezo-composite presented 
by Marat-Mendes had very low cost compared with other sensors [17]. A novel micro-fluidic channel 
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angular accelerometer was proposed by Wolfaardt, and the sensor with a resolution of 15 mu rad/s 
and a bandwidth of 50 Hz was demonstrated [18]. 

In this paper, an angular acceleration sensor without the range limitation of the rotation angle is 
proposed to measure the instantaneous angular acceleration of the rotating system directly, based on 
the principle of electromagnetic induction. In addition, the installation of the sensor is very convenient, 
as the sensor just needs to be coaxially connected with the rotating system when it is used. The output 
signal of the sensor can be led directly without the slip ring to avoid the weakening of the signal. This 
paper describes the mechanical structure and working principle of this angular acceleration sensor, and 
constructs the mathematical model of the sensor, including the transfer function and state-space model. 
The asymptotical stability of the sensor is verified by the Lyapunov theorem when the sensor is 
working. Finally, the sensor calibration is demonstrated, and the probable existing errors and their 
impact on the applications of the sensor are analyzed. 

2. The Structure and Principle of the Sensor 

2.1. The Mechanical Structure of Sensor 

The 3D assembly diagram of the angular acceleration sensor based on electromagnetic induction is 
shown in Figure 1, which includes front cover (T), bearing (2), sensor shaft (3), cup-shaped rotor (3), 
case (5), outer stator core (6), output windings (7), excitation stator core (8), excitation windings (9), as 
well as the end cover @. The excitation stator core (8) is composed of a high magnetic permeability 
soft magnetic iron-nickel alloy sheet or high magnetic permeability silicon steel through punching and 
shearing, and the excitation windings (9) are embedded in its grooves; excitation windings (9) are used 
for the generation of the constant air-gap magnetic field; the cup-shaped rotor (3) is a thin-walled 
non-magnetic cup, which is prepared from high resistivity phosphor bronze, silicon manganese bronze 
and tin zinc bronze and other materials; the sensor shaft (3) is connected coaxially with the shaft of the 
measured rotating system when measuring angular acceleration, and cutting the air-gap magnetic field 
formed by the excitation windings (9); the outer stator core (6) is made of the same materials as the 
excitation stator core (8), and output windings (7) are embedded in its grooves; the output windings (7) 
will generate the electrical signal which corresponds to the angular acceleration of the rotating testing 
system. The corresponding 3D sectional view of the sensor is shown in Figure 2, and the sectional 
sketch map of the sensor is shown in Figure 3. 



Figure 1. 3D assembly diagram of the angular acceleration sensor. 
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Figure 2. 3D sectional view of the sensor. 





Figure 3. Sectional sketch view of the sensor. 




Windings 



2.2. The Principle of Sensors 



The magnetic circuit is shown in Figure 4. When the excitation windings have DC current I\ 
assuming the magnetic circuit is symmetrical, the magnetic flux Oi of the magnetic circuit can be 
expressed by the following equation according to the magnetic circuit theorem [19]: 



2W X I X 



- + - 



-+ 



1=1 



- + - 



l - ) 



= K X I X 



(1) 



In which: 



W\ The effective windings number of the excitation windings 
6\ The effective length through the air-gap <J?i 
S\ The effective area through the air-gap cpi 
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juo The magnetic permeability of the air 

Ir The effective thickness of the cup-shaped rotor 

fiR The magnetic permeability of the cup-shaped rotor 

Sri The effective area through the cup-shaped rotor ®i 

Iji The each segment effective length through the outer stator core of magnetic flux ®i 

Sfi The each segment effective area through the outer stator core of magnetic flux ®i 

fif The magnetic permeability of the outer stator core 

l ti The each segment effective length through the inner stator core of magnetic flux ®i 

S t i The each segment effective area through the inner stator core of magnetic flux ®i 

fii The magnetic permeability of the inner stator core 



- + - 



- + 



+ - 



/u 0 S l MrSrx i=\ M/Sfi Mt^ti 



Figure 4. The magnetic circuit of ®i. 




The cup-shaped rotor can be seen as a number of squirrel-cage rotors. When the magnetic flux is cut 
counterclockwise by the cup-shaped rotor with a speed of n (unit is r/min), the generated motional 
electromotive force is: 



e R =C e Q> 1 n 



(2) 



tjD 



where C e is the structure constant associated with the cup-shaped rotor, C e = , where r is the pole 

moment of the stator and D is the outer diameter of the cup-shaped rotor. The resistance of the 
cup-shaped rotor is r R , and the reactance of cup-shaped rotor can be ignored, the current i R in the 
cup-shaped rotor is: 



i =^ 

L R 



(3) 
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i R generates a magnetic field, and its magnetic circuit is shown in Figure 5. 

Figure 5. The magnetic circuit of O r . 




Similarly, the magnetic flux O r generated by i R can be expressed by the following equation 
according to the magnetic circuit theorem: 



<■>* = 



W R i R 



— + — 1 — + 



E(-V+^-) 

i=i Mf^n Mt^n 



(4) 



In which: 

Wr 

S R 

Sr 

S R 2 

hi 

Sfi 

In 
Sn 



The equivalent number of windings for the cup-shaped rotor 
The effective length through the air-gap 
The effective area through the air-gap 

The effective area through the cup-shaped rotor of the magnetic flux 
The each segment effective length of magnetic flux O r through the outer stator core 
The each segment effective area of magnetic flux O r through the outer stator core 
The each segment effective length of magnetic flux O r through the stator core 
The each segment effective area of magnetic flux through the stator core 
W 9 



8, 



- + - 
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z 
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Fi 



- + - 



fJ f S Fi /d t S Ti 



The total effective windings of the output windings are W2. According to the principle of 
electromagnetic induction, the induction electromotive force in the output windings induced by the 
magnetic flux O r is exhibited as the following equation: 

d(2Q> R ) 



dt 



(5) 



By the combination of Equations (1-5), the following equation can be obtained: 
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-2K x K R C e W 2 W R I x dn 
r B dt 



(6) 



From the above equation, we can find that the induction electromotive force of output windings e Q 
is proportional to the rotation angular acceleration of the cup-shaped rotor. 

3. The Mathematical Model of Sensor 

3.1. Transfer Function 

The equivalent circuit of the working sensor is shown in Figure 6: 

Figure 6. The equivalent circuit of the working sensor. 
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For the cup-shaped rotor, the voltage balance equation is: 

di D 



e R ^R r R + L R 



dt 



For the output windings, the voltage balance equation is: 

e Q =i 2 r 2 +L — + i 2 Z 
at 



(7) 



(8) 
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The output voltage of sensor, namely the voltage across load is: 



U 2 — 



(9) 



The Laplace transform is demonstrated for the Equations (5-9), after combination, the transfer 
function of sensor is obtained: 

L[u 2 ] _ u 2 (s) _ -C e K x K R W 2 ZI x 

(10) 



G(s) 



L [™ ]) sn(s) (r R + L R s)(r 2 + Z + L x s) 
dt 



Jr (r +Z) 
and the 



L R L X 



(r + Z)L +r L 

damping ratio is ^ = , 2 = x for the working sensor system, respectively. The poles of the 

2^{r 2 +Z)L R -r R L x 

T V H~ Z 

system are s x = — ^and s 2 = — respectively. Owing to the fact that the load resistance is large 



L 



and the output windings leakage inductance L x is small, the dominant apex of the system is si, that is 
2C K K W ZI 

G(s) « — — ^—^ — - . The system can be approximated as a series of the proportional components 

r R +L R s 

and the inertial element, and the time constant of the inertial element is . High resistivity material 
can be used to prepare the cup-shaped rotor, in order to improve the response speed of the system. 

3.2. State-Space Model 

Appropriate state variables are selected, the diagonal form of the system state-space model can be 
obtained through the transfer function of the system (10): 
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A. 
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The corresponding state variables scheme is shown in Figure 7: 

Figure 7. The state variables scheme of sensor. 
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The stability of the system can be divided into external stability and internal stability. The external 
stability is defined as the influence of the external input on the response of the system at the initial 
state, which can be judged according to the equation root situation of the closed-loop characteristic 
equation for the system. The internal stability is only related to the impact of the system itself at the 
initial state, also called Lyapunov asymptotical stability. 
Pi Pi 



Assuming P = 



Pi Pa 



the state matrix is introduced in Lyapunov equation, the solution is: 



Pi 



2r D 



P 2 =P3= 0 

Pa 



(12) 



2(r 2+ Z) 



Owing to p\ > 0 and 



Pi Pi 
P3 Pa 



4(r 2 +Z)r R 

sensor is Lyapunov asymptotic stable over a wide range. 



> 0 , the matrix P is positive definite, so the working 



4. The Sensor Calibration 



4.1. Composition of the Calibration Equipment 



The calibration equipment of the sensor is shown in Figure 8, which is composed of bracket CD, 
angular acceleration sensor (2), couplings (3), mass @, coil spring (5), torsion bar (6), base (7), angle 
sensor (8), digital oscilloscope (9), and DC power @. 



Figure 8. The calibration equipment of angular acceleration sensor. 
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The bracket is used to fix angular acceleration sensor (2) and high-precision angle sensor (8); The 
coupling (3) is used for the coaxial connection of the angular acceleration sensor (2) and high-precision 
angle sensor (8); The mass (3) is a self-made iron hoop with holes and has variable quality for selection 
and combination. The moment of inertia for the calibration system can be changed to realize 
calibration at different frequencies through selecting and combining masses; the combination of a large 
mass and a small mass is shown in the Figure; the coil spring (5), the torsion bar (6) and the base (7) 
form a torsional pendulum, which is the power source of torsional vibration for the entire system. The 
torsion bar (6) passes the mass and connects with the angular acceleration sensor ® through the upper 
end of the coupling (3), and its bottom end is connected with the end of the angle sensor (8); The dual 
channel input of the digital oscilloscope (9) connects with the angular acceleration sensor and the angle 
sensor, respectively, for the acquisition and preservation of the output waveform of the two sensors 
during the experiments. 

4.2. The Principle of the Calibration Equipment 

An initial force is applied to the torsion bar @, resulting in the coil spring (5) rotating an angle 6 to 

produce the restoring moment T\ s = kO, in which k is the stiffness coefficient of the spring; The system 

begins to rotate under this restoring moment when the initial force disappears. In the rotation process, 

d6 

the system is affected by the air resistance torque T ls = —y — , in which y is the drag torque 

dt 

dO 

coefficient; and also affected by the bearing friction resistance torque T 3s = -a — , in which a is the 

dt 

friction torque coefficient; Assuming the system moment of inertia is J, we can obtain the following 
equation according to the torque balance equation: 

d 2 6 d6 dO 

— Y + a — + y — 

dt dt dt 



J^Y + a — + y — + kO = 0 (13) 



iy _|_ 

Assuming 2p = - and co 2 n = — , where /? is the damping coefficient of the system, co n is the 

J J 

intrinsic oscillation angular frequency of the system: 



d 2 0 „„d6 



Aj+2 ,_ +atffl = 0 (,4) 

Q 

The solution for above equation is 9{t) = — r =^^=e ~ pt cos(co d t-<p) , in which 0 O is the initial 

amplitude of the pivot angle, co d =co n ^\-^ 2 is the damped oscillation angular frequency, damping 
y B 

ratio is £ = — -j= = — . It is found that the vibration amplitude decays with time; the bigger the 

2-4 kJ co n 

damping factor is, the faster it will decay. Assuming T is the cycle damping vibration for the 

Q 

pendulum, — = e pNT can be obtained if the initial amplitude Oo and the amplitude 9 N of the N-th cycle, 

On 

as well as the used time of N cycles are measured. The damping coefficient can be expressed as: 
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/? = ^-hA 
NT 0„ 



In 

The rotation cycle of the system is T = — = 



2;r 



In 



k An 2 o2 
then — = — ' °~ 

2 J T 2 



(15) 

+ p l can be 



dO 



inferred. The angular rate of the pendulum is — = 0 when system at the amplitude position. At this 

dt 

time, Equation (13) can be simplified as: 



d 2 e _ k e _ 

dt 2 ~ J 



(16) 



where /? is calculated by Equation (15) through experiments, and then the cycle T and the rotor angle 9 
are obtained through the output waveform of angle sensor saved by the oscilloscope. From 
Equation (16), we can get the actual angular acceleration value of the system, corresponding to the 
output voltage of the angular acceleration sensor. Additionally, we can also obtain the sensitivity 
coefficient of the angular acceleration sensor. 

4.3. The Calibration Experiment Results 

The DC resistance of the excitation windings for the angular acceleration sensor is 838 Q. After a 
long-time test at the working temperature, excitation windings voltage U\ = 40 V is selected for the 
angular acceleration sensor. In this case, the temperature of the sensor is about 25 to 28 Celsius 
degrees; the initial angle is about 180° when the external force is applied on the coil spring, the output 
signals of angular acceleration sensor and the angle sensor are shown in Figure 9 and Figure 10 after 
the external force disappears. The output electrical signal of the angular acceleration sensor is 
de-noised by a Rigrsure wavelet; the obtained signal is shown in the Figure 1 1 . The calculated results 
are shown in Table 1 by selecting the signal peaks of angular acceleration sensor and the angle sensor, 
according to the relevant formulas. 

Figure 9. Output signals of angular acceleration sensor. 
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Figure 10. Output signals of angle sensor. 
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Figure 11. The de-noised signals by Rigrsure wavelet. 
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Matlab is used to fit the data in Table 1 and the obtained results are shown in Figure 12. It is found 
that the sensitivity coefficient of the angular acceleration sensor is about 17.29 mv/(Krad/s 2 ). 



Table 1. Experimental calibration data. 



The peak voltage of the angular 


7.32 


7.15 


6.44 


5.44 


5.06 


4.48 


4.04 


3.61 


3.09 


acceleration sensor (mV) 




















The peak voltage 
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of angle sensor (V) 




















The time corresponding peak (s) 
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3.656 
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Oscillation period T/(s) 
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+ 0.494 + 
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Angle values (rad) 
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System damping ratio coefficient (3 
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Figure 12. Calibration curve of angular acceleration sensor. 



7.5 r 

> 

« 7 




?50 200 250 300 350 400 450 

Angular Acceleration (rad/s2) 



4.4. Analysis of the Calibration Experiments 

There will exist a phase difference when the angular acceleration sensor and the angle sensor are 
working in practice. From the calibration principle, we can see that the peak value of the angular 
acceleration is the corresponding amplitude position of the pendulum. Therefore, the peak voltages of 
the two sensors are needed in the calibration process and the phase difference shows little impact on 
the calibration results. Additionally, the aerostatic bearing can be considered to reduce the influence of 
bearing friction on calibration results; the calibration range of the system is limited and only suitable 
for the calibration of the angular acceleration in the lower frequency owing to the constraint of the 
mechanical structure and the limitation of variation range, even though the calibration frequency can 
be changed by the mass to a certain extent; the obtained calibration frequency is not continuous 
because the calibration system depends on the mass. 

In order to solve the issues of non-continuous calibration point and lower frequency band, an 
auxiliary power source can be added to the electromagnetic force to realize a forced vibration state 
when the calibration equipment is working. The linear regulator of the vibration frequency can be 
realized by changing the rotational speed of the power source, resulting in the extension of the 
calibration frequency. These improvement ideas will be demonstrated in our future experiments. 

5. Error Analysis 

5.1. Fluctuations of Excitation Voltage 

From the sensor principle, it is recognized that magnetic flux Oi generated by the excitation current 
includes the variation component AOi if the excitation voltage has fluctuations. The cup- shaped rotor 
cuts the magnetic flux Oi with a constant rotation speed, so its motional electromotive force eR 
includes the variation component Ae R , resulting in the existence of a current change component Ai R for 
the cup-shaped rotor, and the generated magnetic flux O r also contains a variation component AO R . At 
this time, output windings will generate the corresponding induced electromotive force, known as the 
residual voltage. In order to overcome the abovementioned disadvantages of the sensor, a constant 
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current source or a permanent magnet can be used to replace the existing excitation circuit to guarantee 
the constancy of the magnetic flux ®i. 

5.2. Effect of Temperature 

The properties of the excitation windings, the output windings, the resistance of the cup-shaped 
rotor and the magnetic material may change because of the change in ambient temperature and the 
heating from extended use of the sensor. These property changes show influences on the output 
induced electromotive force of the sensor, resulting in the instability of the output performance. The 
character of the temperature-sensitivity coefficient of the sensor is shown in Figure 13. In this 
experiment, DS18B20 is used to detect the temperature of the sensor case. According to the results, the 
sensor sensitivity coefficient will decrease slightly with increasing temperature. 

Figure 13. Temperature-sensitivity coefficient characteristics of the sensor. 

» 18.5, , , , , , , , , , 1 

to 




^ 16 10 15 20 25 30 35 40 45 50 55 60 

Temperature / a C 

In order to decrease the impact of temperature, temperature compensation should be applied. A 
negative temperature coefficient thermistor can be concatenated in the excitation and the output circuit 
to compensate the temperature variation influence. 

5.3. Nonlinearity Errors 

The cup-shaped rotor of actual sensor is a first-order circuit composed of a resistance and 
inductance, so the change of the current ir lags induces an electromotive force eR due to the influence 
of the inductance Lr, resulting in the distortion of the symmetrical magnetic circuit: the magnetic flux 
Or generated by the current i R varies with the change of rotational speed, and the change component of 
Or demagnetizes magnetic flux ®i, causing the change of the magnetic flux ®i. As a result, linear 
relationships of the output potential and the angular acceleration are destroyed, resulting in the 
presence of a non-linear error. 
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6. Conclusions 

The angular acceleration sensor presented in this paper shows strong practicality without limitation 
on the rotation angle range. It only needs to connect with the measured rotating system coaxially when 
working. The structure and principle of the calibration equipment are described, and the calibration 
curves and performance indicators of the sensor are obtained by experiments. In our following work 
efforts will be made to examine the practical applications for this sensor in view of the extensive 
applications of angular acceleration sensors, such as transmission noise measurement and fault 
diagnosis for complex gearboxes or the internal combustion engine, fuel injection timing and fuel 
injection amount control for automotive engines, as well as the kinetic analysis for mechanical arm 
joints and so on. 
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